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ALGEBRAIC PROPERTIES OF SELF-ADJOINT SYSTEMS* 



DUNHAM JACKSON 

The general definition of adjoint systems of boundary conditions asso- 
ciated with ordinary linear differential equations was given by Birkhoff.f 
In a paper of Bocher, J in which the idea is further developed, there is obtained 
a condition that a system of the second order be self-adjoint. It is proposed 
here to extend the discussion of this problem to the case of systems of any 
order. § A condition for self-adjointness of the boundary conditions is simply 
expressed in matrix form, without any requirement that a corresponding prop- 
erty be possessed by the differential equations; the condition gains in sym- 
metry if it is assumed that the differential expression which forms the left- 
hand member of the given differential equation is itself identical with its 
adjoint or with the negative of its adjoint. 

As a preliminary, it will be well to recall a well-known rule for the combina- 
tion of matrices, which will be found particularly convenient. Let a\, a 2 , 
az , and a 4 be re-rowed square matrices. The square matrix of order 2n which 
is made up of the elements of these four, with the elements of a\, arranged 
in order, in its upper left-hand corner, and the elements of the other matrices 
correspondingly disposed, may be indicated by the notation 



\ «3 OLi ) 



* Presented to the Society, December 27, 1915. 

t G. D. Birkhoff, these Transactions, vol. 9 (1908), pp. 373-395; p. 375. 

JM. Bdcher, these Transactions, vol. 14 (1913), pp. 403-420; p. 409. Cf. E. 
Bounitzky, Journal de math^matiques pures et appliqu^es, 6th 
series, vol. 5 (1909), pp. 65-125; pp. 112-116; W. Stekloff, Sur certaines questions d'analyse 
qui se rattachent a plusieurs problemes de la physique mathematique, Mgmoires de 1 ' A - 
cademie Imp6riale des Sciences de St.-Pltersbourg, classe physico- 
math^matique, 8th series, vol. 31 (1913), no. 7; pp. 5-6, 69. 

§ Bounitzky, in the paper just cited, pp. 107-116, gives a discussion of the problem, laying 
stress primarily on the case of a self -adjoint differential equation. His method of treatment 
is different from that employed here, and in fact bears as little resemblance to it as could 
well be expected, considering that the author's task in either case, as far as this particular 
problem is concerned, is merely to give a concise account of algebraic processes which in them- 
selves are perfectly straightforward. 
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Let /3 be another matrix of the same sort, consisting of the elements of the 
n-rowed matrices ft , • • • , (84 . The notation suggests the symbolic product 

fa x a 2 \/ft ft \ = / on ft + « 2 ft aift + a 2 ft\ 
\a 3 a4/\ft (84/ \a 3 ft + a4ft a 3 ft + a4ft/' 

formed as if the elements indicated were ordinary numbers instead of matrices. 
By actually writing out the n-rowed matrices «i ft + « 2 ft , etc., it is possible 
to interpret the product as a square matrix of order 2n; and it is found that 
this is the same as the product obtained by multiplying together a and j3 as 
2n-rowed matrices in the usual way. 

To begin the discussion of the problem in hand, let L(u) stand for the 
differential expression 

T d n u . . d n_1 u . 

L(u) = MaO-^ + Pi(aO-^=r+ ••• +p n (x)u, 

where n is any positive integer. The coefficients p are real or complex func- 
tions of the real variable x , defined in an interval a Si x Ss b , and it may 
be assumed, with more regard for simplicity of statement than for ultimate 
generality, that they are continuous together with their first n derivatives. 
With the differential equation 

L(u) = 

there is supposed given a set of n linearly independent boundary conditions 
U,(u) = 0, i — 1,2, • • • , n, in which the general left-hand member is a 
linear combination, with real or complex constant coefficients, of the 2n 
quantities u(a) , u' (a) , • • • , w (n-1) (a) , u{b) , u' (b) , •••, u^"^ ( b ) . Let 
n of these quantities be chosen in such a way that the determinant of the 
corresponding coefficients is different from zero. The possibility of such a 
choice is a consequence of the assumption that the n conditions are linearly 
independent. The quantities selected will be denoted by Ui , u% , • • • , u n , 
the remaining n quantities by u„+i , • • • , ui n ; the order of the subscripts within 
each group is immaterial. It is most natural to assign the first n subscripts 
toit(o), • • • , u < - n ~ t> (o), or to «(&), •••, w (n_1) ( b ) , if the corresponding 
sets of terms are linearly independent, but this will not always be the case. 
The differential expression adjoint to L (u) is 

M{V) = (-l)n^ n{poV ) + (-l)»-l ^- i (p lV ) + ... +p nV . 

It is connected with L ( u ) by Lagrange's identity, 

vL(u) — uM (v) = -j- P (u, v) , 
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where P(u,v) is a bilinear form in u(x) , u' (x) , •••, « (re-1) (a;); v(x), 
»' ( x ) , • • • , » (re_1) ( x ) , with coefficients which are functions of x . Integrated 
with regard to x from a to b, this identity gives a relation involving the values 
which the functions just named take on at the ends of the interval. If the 
values of v and its derivatives at a and b are denoted by Vi , • • • , v 2n , in such 
a way that Uj and vj refer to derivatives of the same order taken at the same 
point, the relation becomes 



f 

•J a 



\vL{u) — uM (»)] dx = II («,«), 



in which II ( u , v ) is a non-singular bilinear form in the 4n quantities Uj , Vj , 
with constant coefficients.* 

The definition of the adjoint boundary conditions proceeds now as follows. 
In addition to the linear forms Ui ( u ) , ■ • • , U„(u) , there are defined n other 
forms Un+i ( u ) , • • • , \J in ( u) , which are arbitrary, except for the restriction 
that the whole set of 2n forms shall be linearly independent. Then In linear 
forms V\ ( v ) , • ■ ■ , V« n ( v ) , in the variables v } - , are determined so that 

In 

n(tt,«) = Z Ui(u)Vi(v). 

The last n of these forms, set equal to zero, express the adjoint conditions. 
It is well known that a different choice of the arbitrary forms U n +i , ■ ■ • , U 2n , 
would lead to a set of conditions equivalent collectively to those already 
obtained, that is, a set of independent linear combinations of them.f It is 
to be kept in mind throughout that it is the set of equations that is important, 
and not the individual equation, still less the individual linear form. 

The definition just given can be expressed in slightly different words, as 
follows: The variables u are subjected to a linear transformation u' ( = Ui(u) , 
and then a linear transformation of the other set of variables, v'i = F,-(i>), 
is determined so that the two transformations together carry over the bilinear 
form II ( u , v ) into the normal form ^jUi v'i . Let II denote the matrix of the 
coefficients in the form II ( u , v); the double use of the letter II , to represent 
a matrix and to represent an operator, will not cause any confusion. Let a 
denote the matrix of the coefficients in the inverse transformation, by which 
the u's without accents are expressed in terms of the u's with accents, and 
let j8 be the corresponding matrix for the v's; the conjugate of a will be de- 



* It is to be noticed that the coefficient p„ is without effect on the form of II , as the term 
p n uv cancels from the difference vL (u) — uM ( v ) . Consequently, if the differential equa- 
tion with a parameter, L (u) -\-\u =0, is under discussion, the form II is independent of 
the parameter, and the same is true of the adjoint boundary conditions, which depend on n 
for their definition. 

f Birkhoff, loc. cit., p. 375; Bocher, loc. cit., p. 405. 
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noted by a' , and, in general, conjugate matrices will be indicated by accents. 
Finally, let I denote the unit matrix of order 2n, having unity for each ele- 
ment in the principal diagonal, and zero everywhere else. Then the relation 
connecting these matrices is simply the equation* 

(1) a'TL(3 = I. 

Tt is to be noticed that not a, but a -1 , is to be regarded as given directly, 
and /3 -1 is the matrix that is ultimately desired. An expression for the latter 
is found by multiplying both members of (1) on the right by/3 -1 : 

/T 1 = a' n . 

Hence the problem to be solved may be formulated thus: 

Under what circumstances will the last n rows of a' II be independent linear 
combinations of the first n rows of aT 1 ? 

Let «i stand for the square matrix of the coefficients of the first n variables 
Uj in the n given forms Ui , and a 2 for the matrix of the remaining coefficients 
in these forms. Since the notation has been supposed chosen so that a\ 
is non-singular, the 2n forms U\, • • • , Uin. will be linearly independent if the 
last n of them are set equal to u n +i , • • ■ , w 2n respectively. Then the matrix 
a -1 has the form 



(on <x 2 \ 



here I denotes the unit matrix of order n (the use of the same letter to repre- 
sent unit matrices of different orders will not cause any confusion) and stands 
for a matrix composed entirely of zeros. It may be verified at once that 

/ ai a 2 \ ( aT 1 — aT 1 a 2 \ _ / 1 0\ 

lo ijlo i )-\o IJ' 

from which it follows that the second factor on the left is the matrix a . The 
product af 1 a 2 will be denoted for brevity by the single letter S . It will be 
convenient also to indicate the conjugate of the reciprocal of a non-singular 
matrix by a double accent, so that (a^ 1 )' = (ai) -1 = a" , and 5' = a 2 a" . 
Then we have for the conjugate of a the representation 



' ( a " °\ 



The matrix of the bilinear form will be represented for the present simply by 
the notation 



\TT4 7T 2 / 



' Cf., e. g., Bdcher, Introduction to Higher Algebra, p. 115. 
Trans. Am. Math. Soc 28 
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which is applicable to any square matrix of order 2m; with regard to the special 
form of II , more will be said later. The expression for /3 -1 is 



M-rl 



TTi + TTi — 5 7T3 + 7T2 , 



Now, since a\ is non-singular, it is possible to express any sets of n quantities 
each, in particular, the rows of the matrix standing in the lower left-hand 
corner of a' II, as linear combinations of the rows of a\\ and the coefficients 
by which this is accomplished are uniquely determined. The question is, 
then, whether the corresponding linear combinations of the rows of a^ are 
identical respectively with the rows of the matrix in the lower right-hand 
corner of a' II . It will be seen presently that if this is the case, the combina- 
tions in question will necessarily be independent. 

In general terms, if n is a matrix of order n , and there are formed n linear 
combinations of its rows, the coefficients in these combinations forming the 
successive rows of a matrix y , the resulting sets of n quantities each are 
respectively the rows of the product matrix 777 . In the present instance, 
7 is to be determined so that 

(2) yai = — 8' iri + in , 

and the condition for a self-adjoint system is that the matrix y so defined 
satisfy the further relation 

(3) ya2 = — 5' 7T 3 + 7T2 . 

If it does, then 

_ / a" 7Ti a" 7T 3 \ _ / I \ / a" in d[ t 2 \ 
V yai Y«2 / \0 )/\ «i «2 / ' 

If 7 were singular, the matrix a' II would be singular, which is certainly not 
the case; it follows that the determinant of y is necessarily different from zero. 
Let (2) be solved for y , and the result substituted in (3). It is found that 

(4) 7 = — 5' WlO^ 1 + 7T4 ai" 1 > — 5' TTl 5 + 1T4 5 = — 5' 1T3 + 7T2 . 

In the last equation, 5 has been substituted in the left-hand member for its 
equal, a^ 1 a 2 • The relation (2), defining y , may be regarded as valid in any 
case. Then either of the relations (3), (4), is a consequence of the other. 
The conclusion may be stated as follows : 

Theorem. A necessary and sufficient condition that the given set of boundary 
conditions be self -adjoint, is that* 

(5) 5' 1Ti 5 — 7T4 5 = 5' 1T3 — 7T2 . 



* Professor W. A. Hurwitz points out that this condition can be written concisely in the 
form A' IIA = 0, where 
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In the cases that are of most interest, the form of the condition can be 
simplified. Suppose first that in the n given boundary conditions, the sets 
of terms corresponding to one end of the interval are linearly independent, 
so that the values of u and its derivatives at this point may be represented 
by Mi , • • • , ti„ . Then, as no single term of II ( u , v ) is concerned with both 
ends of the interval, the subscripts in each term that actually appears will 
be both less than n + 1 or both greater than n , and the matrices previously 
denoted by ir 3 and 7r 4 will be composed entirely of zeros.* If 5 is replaced by 
the alternative expression a^ 1 a 2 , the condition for self-adjointness reads as 
follows: 

«2 <xi tti aT 1 «2 = — Ti. 

It appears that in this case ai as well as a.\ must be non-singular, since -wi is 
non-singular. On multiplication before and after by cl' 2 and a^ 1 respectively, 
the condition takes the form 

Oil 7Ti Of' = — Oi'z 7T2 Oi^ 1 , 

the meaning of which can be expressed in still different language. There 
is.no loss of generality in assuming, for definiteness, that on is the matrix of 
the coefficients belonging to the point b . More precisely, it may be assumed 
that 

Mi = u ( b ) , u 2 = u' ( b ) , • • • , u„ = w (n_1) ( b ) , 

w n+1 = w(a), u n +2 = u'(a), •••, u 2n = w (n_1) ( a ) , 

while Vi , • • • , Vi n have a corresponding interpretation. Let the value of 
P(u, v) for x = b be denoted by Hi ( u , v ) , and its value for x = a , by 
n 2 ( u , v ) , so that 

II ( u , v ) = IIi ( u , v ) — n 2 ( u , v ) . 

The expressions IIi and II2 are bilinear forms in u\ , • • • , u n ; «i , • • • , v n , and 
Un+it •••, u in ; fl n +i, •••> «2n, respectively, and their matrices are tv\ and 
— 7T2 . Let each of the n given linear forms Ui(u) be regarded as the sum 
of a form Un ( u ) , involving only the first n of the variables u, , and a form 
Uii{u), involving only the last n variables. If Hi(u,v) is written as a 
bilinear form in Un (u) and Ua (v) , i = 1,2, • ■ • , n, the coefficients in the 
representation have the matrix a,' t 1 af l . A corresponding remark applies 
to II2 ( u , v ) . It follows that 

If the forms Un (u) are linearly independent, the condition for a self-adjoint 
si/stem is that the coefficients by means of which Hi (u, v) is expressed in terms 
of the linear forms Un ( u ) , Un ( » ) , be the same as those by which II2 ( u , v) 
is expressed in terms of the forms Ua ( u ) , Ua ( v ) . 

* The fact that in and V2 also contain a considerable number of zeros need not be taken into 
account. 
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Now let the restriction on the choice of the first n variables be removed, 
but let it be supposed that when the given differential expression L ( u ) and 
its adjoint are written with the same function u as argument, the latter ex- 
pression is identical with L ( u ) or with the negative of L ( u ) , in other words, 
that L(u) is self -adjoint or anti-self -adjoint. The statement of the theorem 
is again capable of some simplification, for the reason now that the bilinear form 
P(u,v) is skew-symmetric in the one case and symmetric in the other,* 
and the same thing consequently is true of the matrix II . More generally, 
it suffices for the symmetry of P and II that L ( u ) differ only in the term p n u 
from an expression which is anti-self-adjoint,f since, as has already been 
remarked, this term has no effect on the form of P . 

Let the condition (5) be written in the equivalent form 

(6) 5' TTi 5 - 2ir 4 S + ir 2 = - 5' tti 5 + 25' tt 3 - tt 2 . 

If II is skew-symmetric, then tc'\ = — wi , r' 2 = — -k% , and v\— — 7r 3 , and 
it appears that the matrix which forms the right-hand member of (6) is the 
conjugate of the matrix on the left. If II is symmetric, each member of (6) 
is the negative of the conjugate of the other. That is: 

7/ the differential expression L(u)is self-adjoint, the condition that the boundary 
conditions be self-adjoint ?'s that the matrix 

5 7!"1 5 — 2lTi 5 -)- 7T2 

be symmetric; if L(u) is an anti-self-adjoint expression, or differs from such 
an expression only in the term of order zero, the condition is that the matrix just 
written down be skew-symmetric. 

This form of statement remains valid, of course, under the earlier hypothesis 
that 7T4 = 0; but in that case the condition simply resumes the form 

5' 7Tl 5 + 7T2 = , 

since this matrix is required to be symmetric and skew-symmetric at the 

same time.J 

Harvard University, 
Cambridge, Mass. 

* Cf. L. Schlesinger, Handbuch der Theorie der linearen Differentialgleichungen, vol. 1, 
pp. 70-71, 73-74. 

f This observation is not trivial, as the corresponding statement with regard to self-adjoint 
expressions would be; for example, the expression xu' + fit is anti-self-adjoint, while the 
expression xu' is not. 

t It is necessarily skew-symmetric if L (u) is self-adjoint, without regard to the boundary 
conditions, and becomes symmetric if the latter are self-adjoint, and similarly in the anti- 
self-adjoint case. 



